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We employ observationally determined intrinsic velocity widths and column densities of 
damped Lyman-alpha systems at high redshift to investigate the distribution of baryons in 
protogalaxies within the context of a standard cold dark matter model. We proceed under 
the assumption that damped Lyman alpha systems represent a population of cold, rotationally 
supported, protogalactic disks and that the abundance of dark matter halos is well approximated 
by a cold dark matter model with critical density and vanishing cosmological constant. Using 
conditional cross sections to observe a damped system with a given velocity width and column 
density, we compare observationally inferred velocity width and column density distributions to 
the corresponding theoretically determined distributions for a variety of disk parameters and CDM 
normalizations. In general, we find that the observations can not be reproduced by the models for 
most disk parameters and CDM normalizations. Whereas the column density distribution favors 
small disks with large neutral gas fraction, the velocity width distribution favors large and thick 
disks with small neutral gas fraction. The possible resolutions of this problem in the context of 
this CDM model may be: (1) an increased contribution of rapidly rotating disks within massive 
dark matter halos to damped Lyman-alpha absorption or (2) the abandoning of simple disk models 
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within this CDM model for damped Lyman-alpha systems at high redshift. Here the first possibility 
may be achieved by supposing that damped Lyman-alpha system formation only occurs in halos 
with fairly large circular velocities and the second possibility may result from a large contribution 
of mergers and double-disks to damped Lyman-alpha absorption at high redshift. 
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1 Introduction 


It is widely believed that the phenomena of damped Lyman-alpha (hereafter; damped Ly-a) 
absorption in the spectra of high redshift quasars is due to large amounts of neutral gas in 
dark matter halos which intervene on the lines of sight to quasars |Wolfe et al. 198^ . It is 
intriguing that the ratio of neutral gas mass density to the critical density in damped Ly-a 
systems, ^^Lyaiz), at redshift 2 ; ~ 3 is comparable to the baryonic content of stars in galaxies 
at the present epoch [Lanzetta et ah 19951 IWolfe et ah 1995j . Both the steady decline in Q.Lya{z) 
with decreasing redshift and the agreement between ^Lya{z) at low z and the neutral gas content 
of nearby galaxies |Rao &: Briggs 1993| are commonly interpreted as the incorporation of baryons 
into stars within the damped systems [Lanzetta et ah 19951 IWolfe et ah 1995] , Observations of 
damped Ly-a systems have, and continue, to yield important information about the abundances 
and masses of dark matter halos, as well as the process of star formation |Ma and Bertschinger 1994| 
IKauffmann and Oharlot 19941 iMo fc Miralda-Escoude 1994] , 

Whereas the hosts of damped Ly-a systems seem to be established, there is considerable debate 
about the distribution of neutral gas in individual dark matter halos at high redshift. For example, 
neutral gas in damped Ly-a systems may be distributed spherically symmetric within the halo, 
or may be within a rotationally supported protogalactic disk. The main observational support for 
damped Ly-a systems representing halo gas comes from a comparison of metallicities of damped Ly- 
a systems as a function of redshift with the metallicities of disk and halo stars as a function of stellar 
age |Pettini et al. 19M1 ILu et al. 1996j . On the other hand, there is accumulating observational 
evidence that damped Ly-a systems at high redshift are protogalactic disks. Prochaska &: Wolfe 
(1997a; hereafter PW) recently analyzed a sample of 17 damped Ly-a systems observed at high- 
resolution with HIRES on the Keck I telescope. They argue that observationally determined velocity 
profiles of high-redshift damped systems are reflective of rapidly rotating, cold disks, and cannot be 
reconciled with isothermal halo gas models or slowly rotating “hot” disks. In addition, the recent 
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observation of the stellar and continuum emission of a high-redshift galaxy believed to be responsible 
for damped Ly-a absorption in the spectrum of a quasar may be explained straightforwardly only 
if one assumes the galaxy to be host of a disk [Djorgovski et al. 1996|ll7u et al. 1997] . 


Unfortunately it is not possible to easily obtain a resolution to the debate about the spatial and 
velocity distribution of neutral gas in dark matter halos from numerical simulations. This is partially 
due to limited resolution, such that simulations which contain enough halos to be statistically 
representative may not resolve the physics within an individual halo [Qarduer et al. 1997] . 
Furthermore, when simulating the formation of disks in individual halos within the context of 
currently popular structure formation scenarios, the sizes of disks are predicted to be much 
smaller than the observationally inferred sizes of present day galaxies, casting significant doubt on 
the predictive power of the simulations Navarro, Frenk &: White 1995| . Nevertheless, numerical 
simulations may shed further light on the nature of damped Ly-a systems; for example, they may 
test if it is realistic to employ simple models for the velocity distribution of gas or if, rather, 
the occurrence of damped systems results from cold, rotating disks, as well as halo gas, merging 
galaxies, and/or irregular galaxies in almost equal parts IFlaehuelt et al. 1997] . 


Semianalytic models for the formation of disk galaxies are very successful in 
explaining properties of present-day spiral galaxies [Fall fc Efstathiou 19801 IKauftmann 19961 
IDalcanton et al. 19971 IMo et al. 1997] . These models consist of two parts, a Press-Schechter type 
formalism to determine the abundances and masses of halos in a given hierarchical structure 
formation scenario, and a model for the formation of a rotationally supported disk within the dark 
halo. The sizes and profiles of disks may be estimated by the distribution of halo angular momentum 
under the assumption that there is no angular momentum transport between halo and disk, and 
that baryons are initially in solid body rotation. When combined with efficiencies to turn gas into 
stars, such models successfully explain the observed trends and scatter in exponential disk scale 
length, surface brightness, and the Tully-Fisher relationship [Dalcanton et al. 1997IIMo et al. 1997j . 
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It is instructive to see if such disk models may also reproduce the observational data of damped 
Ly-a systems at high redshift. In fact, an initial, not fully self-consistent, comparison between the 
observational data and one particular disk formation model |Ka,nffma.nn 1996] demonstrated that 
there is significant discrepancy between the model and the observations (PW). 

In this paper we compare observationally inferred velocity widths of a sample of damped Ly-a 
systems at mean redshift 2 : ~ 3 to theoretically predicted velocity widths for a variety of disk 
models. In Section 2 we summarize the assumptions entering our calculation of halo abundances 
and masses, as well as disk rotational velocities. We also describe our computation of conditional 
cross sections, to observe a disk with more than the critical column density for a damped system 
and a velocity width more than a certain fraction of the rotational velocity of the disk. Our results 
and conclusions are presented in Section 3. 

2 Semianalytic Model of Protogalactic Disk Populations 

We consider the conditional optical depth per unit redshift, dr/dz{Avth), for observing a damped 
Ly-a system showing a velocity spread, Av, larger than some threshold, Avth- Note that the optical 
depth per unit redshift is frequently called the rate of incidence. The conditional optical depth is the 
product of proper number density of halos hosting damped systems, the conditional cross section 
for observing a damped system with velocity spread larger than Avth, and the proper path length 
traveled by a photon per unit redshift interval, i.e. dr/dz = n^addF/dz. We assume that the 
kinematics and gas distribution of damped Ly-a systems is adequately described by rotationally 
supported disks within dark matter halos. To account for a continuous distribution of halo masses 
we use the Press-Schechter formalism for halo formation with density threshold 5c = 1.68 and 
top hat filter. In our calculations we assume that the abundances of halos are well approximated 
by the abundances in a closed cold dark matter model with vanishing cosmological constant and 
varying normalizations. We use the CDM transfer function by Bardeen et al. (1986) with scale 
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invariant initial perturbations (n = 1) and Hubble parameter Hq = 50 km s ^ Mpc ^ {h = 0.5). 
The conditional optical depth for a continuous mass distribution may be written as 


= ^(1 + r Au*;,), 

dZ tlQ J (LlVl 


( 1 ) 


with c the speed of light. Note that we will use an equation similar to Eq. ^for the computation 
of column density distributions where the conditional disk cross section in Eq. ^ is replaced by 
an adequate column density distribution. In writing Eq. ^ we assume that damped systems only 
form within sufficiently large halos with circular velocities exceeding a minimum value, u™*”. Even 
though it has been argued on physical grounds that u™*” ~ 35 — 50 km s“^ due to photoionization by 
the UV-background at high redshift and the associated long cooling time scales for baryons in small 
halos [Efstathiou 1992| . we will treat as a free parameter in order to assess the contribution 
of less massive halos to the conditional optical depth. 


We may parametrize the conditional disk cross section by 


(Jd{M, Avth) = TtRI F{Nj_, k\N > Nd,fv = Avjvrot > AvthjVrot) , ( 2 ) 

where Rd, N±, k, and Vrot are scale length, central column density, ratio of scale height to scale 
length, and rotational velocity of the disk, respectively. The cross section receives only contributions 
from lines of sight which (1) result in column densities N exceeding the threshold for a damped 
system Nd = 2 x 10^*^ cm“^ and (2) show velocity width Av = fvVrot exceeding the velocity width 
threshold, Avth- Note that due to the self-similarity of disks with different Rd but the same 
and K, the function F depends explicitly only on and k, but not Rd, provided one assumes a 
flat rotation curve. We assume the disk rotational velocity equals the circular velocity of the halo 
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(3) 

Here r^r is the virial radius of a halo calculated within the spherical collapse model, and G is the 
gravitational constant. Motivated by the successes of semianalytic formation models of present-day 
disk galaxies, we assume that the scale length of the disk is proportional to the virial radius of the 
halo 


Vrot — Vc — 


^ GM km , 1 /n 

-= 159.6- 


1012 Mq 




Rd — vir — 


(179)1/3 \47rpo/ (1 +-s) 


= 169.7kpc /- 


h-2/S / M 
(1 + z) V1012 Mo 


In order to evaluate conditional disk cross sections, we make several simplifying assumptions 
regarding the physical characteristics of the disks. We adopt a simple, phenomenologically 
motivated disk model, the exponential disk. As in PW, we parameterize the neutral gas density of 


the disk, n//, by 


nniR, Z) = riQ exp 




(5) 


where no is the central gas density. Furthermore, we describe the gas kinematics by a flat rotation 
curve, parameterized by the rotation speed, Vrot, and a random motion component characterized 
by a 1-dimensional velocity dispersion, cr = 10 km s “i . Conserving baryonic mass within halos, 
the central column density may be related to the total halo mass by 


N±= [ nniO, Z) dZ = 2noKRd = ^ ( 6 ) 

J nmH2TrR^ 

where mn is the hydrogen mass, = 0.05 is the assumed fractional contribution of baryons to 
the critical density, and ^ ~ 1.3 accounts for the fact that a fraction of all baryons are in form of 
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helium. The parameter /at < 1 is the fraction of baryons participating as neutral gas in the disk. 
This fraction may be smaller than one if either there is an ionized component of gas due to partial 
photoionization or if gas is efficiently incorporated into stars. 

In order to calculate the cross section ad{M,Avth) we consider two conditional cross- 
sections: (1) (Jnhi, defined as the cross-section for sightlines yielding column densities A^(HI) > 
2 X 10^*^ cm“^ (the damped Lya threshold), and (2) the conditional velocity width cross-section, 

for sightlines which give a velocity width, Av, greater then a fraction, fc, of the rotation speed 
(/c e [0,1]). We calculate and a^Hi by systematically running sightlines through an exponential 
disk according to the techniques outlined in §3.2 of PW. While computer intensive, the process is 
straightforward and reproducible. To verify the accuracy of our procedures we compare (Tnhi 
with the cross-sections determined by Fall & Pei (1993) in the thin-disk limit and find excellent 
agreement over the whole column density range = 10^®'® — 10^^'®cm“^. Note also that we 
find the thin-disk limit to be an excellent approximation to ctnhi even for finite thickness disks 
(k ~ 0.5). To simplify the computations, we make use of the fact that is independent of N^. 
We may then obtain the doubly conditional cross section ad for given k and fc by convolving 
the cross-sectional area leading to a^Hi with the cross-sectional area leading to aj^. 

The doubly conditional cross-section depends on the disk parameters in the following ways. 
Thicker disks (larger k values) imply larger velocity widths and therefore larger ad for a given fc 
value due to the increased path length of sightlines within the disk. The dependence of a^ on iV_i_ 
is more complicated. As noted in PW (e.g. Figure 3), the velocity widths tend to decrease with 
larger impact parameters for a given sightline. This is primarily because the differential rotation 
projected along the line of sight decreases with increasing impact parameter. For a disk with a 
large N± value even sightlines at a number of scale lengths may lead to the observation of a damped 
system. However, due to the large impact parameter such sightlines only yield fairly small Av/vrot- 
These ideas are illustrated in Figure 1 where we plot a^, for two different N± = 10^®'®, 10^^'^cm“^ 



and K = 0.108,0.232 as a function of fc, and where ad = ajsrHi for fc < 0.05. It is essential for the 
conclusions of this paper that we find only a small fraction ~ 1/10 of sightlines which result into 
the observation of a damped system also show more than ~50% of the rotation velocity for large 
A^_l disks. This is not the case for small A^_l disks where this ratio is approximately one-half. 

3 Results 


Given the results from §2, we can calculate the optical depth per unit redshift, dr/dz , for damped 
Lya systems showing more than a certain velocity width, Avth, as a function of Avth- The model 
has 5 free parameters: (1) erg, the linear rms fluctuation amplitude measured at 8 h“^ Mpe, (2) 
the circular velocity cutoff corresponding to the least massive dark matter halo capable of 
forming a damped system, (3) k, the ratio of the scale height to the scale length, (4) /, the scale 
length of the disks as a fraction of the virial radius, and (5) /at, the fraction of baryons within a 
dark matter halo which contribute to the disk. In Figures 2a-e we compare optical depth curves 
calculated over a wide range of the physically allowed parameter space with dr/dz measured from 
26 damped Lya systems with redshifts z = 2.0 — 4.4 [z = 2.86). Fourteen of the systems were 
presented in PW and the remaining 12 can be found in Prochaska & Wolfe (1997b). The sample is 
kinematically unbiased and each system has an accurately determined velocity width. 


Figures 2a and 2b present the effects of erg and respectively, on the predicted optical 

depth. These two parameters determine the number density of halos as a function of circular 
velocity. In particular, by increasing cjg and/or t;™"' the mass function of halos capable of forming 
damped systems is shifted towards halos with larger Vc, and disks with larger Vrot- In panel (a) 
we consider a range of erg values, including the two most favored observational values: cjg = 0.667 
from the observed masses of rich clusters [White, Etstathiou, &: Frenk 1993 and erg = 1.2, the 
normalization implied by measurements of the CMBR [Gorski et al. 19^ . First note that none of 
the dr/dz curves in panel (a) are consistent with the observational data. Whereas the observations 
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indicate that a fraction of one-half of all damped systems show an intrinsic velocity width of 
Av > 100 km s“^ our models show that this fraction is < 1/10 for ug < 1. Note that we choose 
/ = 0.1 and /at = 1 in panel (a) yielding typical central column densities of A^_l 10^^ cm The 
small fraction of large Av systems in our model is then due almost in equal parts to a comparatively 
large fraction of small Vc systems as well as due to the small chance for observing a significant 
fraction of the rotational velocity at such high Nj_ (cf. Fig. 1). It is evident from panel (b) that 
by removing the contributions of small halos {vc < 200km s ) to the optical depth one may 
obtain approximate agreement between observations and the model. This highlights the fact that 
the CDM model with ug ~ 1 predicts enough massive halos {vc > 200 —300km s“^) at high redshift 
to account for the abundance of damped systems with large intrinsic velocity width. 

The dependence of the optical depth on the parameters k, f and /w which affect the cross- 
section of the disks within the dark matter halos, is shown in panels (c) —(e). In Figure 2c we 
plot dr/dz against Av for models with varying disk thickness. It is seen that the observational 
data favors thick disks and that disks with thicknesses comparable to the thin disk of the Milky 
Way (k ~ 0.05; Mihalas &: Binney 1981) are unacceptable. We investigate the effect of varying the 
disk size in panel (d). One may produce a nearly acceptable fit of the velocity distribution data 
even for low normalization CDM models (erg = 0.5) by taking the scale length of the disks to be 
a considerable fraction of the virial radius (/ ~ 0.8). Agreement between the model predictions 
for the velocity width distribution and the observations is then attained due to the comparatively 
large typical fraction of the rotational velocity one observes in these low central column density 

r\j 1020.3 

cm ^) systems (cf. Fig. 1). It is important to stress that analytic models of disks 
within a CDM cosmogony (e.g. Dalcanton et ah 1997, Mo et al. 1997) suggest / ~ A, the average 
spin angular momentum parameter which has been measured numerically to be A ~ 0.05 ± 0.03 
[Barnes fc Efstathiou 19^ . Such analytic models are thus inconsistent with extremely large scale 
length Rii- Finally, we allow for the possibility that not all of the baryons within the halos will 
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collapse to form a gaseous disk. Numerical simulations suggest |Ma et al. 1997j that a fraction of 
baryons in halos will remain hot due to photoionization by the UV-background and therefore not 
contribute to damped Lya absorption. Moreover, it is possible that a significant fraction of baryons 
may be incorporated into stars by a redshift of z ~ 2.5. As is evident from panel (e) in Fig. 2, 
lowering /at has the effect of increasing the fraction of high Av systems, thereby slightly lessening 
the discrepancy between model and observations. This is again due to a decrease in the average 
central column densities of disks. 

Nevertheless, disks with low central column densities are inconsistent with the column density 
distribution for damped Lya systems at high redshift. In Fig. 3 we compare the column density 
distribution resulting from several of the models in Fig. 2 to the observationally determined column 
density distribution at mean redshift z = 3.0 [Wolfe et al. 199^ . It is evident from the figure that 
the observational data clearly favors large N± disks (the dotted and long-dashed lines are models 
with low average A^_l). Significantly increased scale lengths Rd, or decreased neutral fractions fjy, 
are therefore ruled out as an explanation for the comparatively flat velocity width distribution of 
damped Lya systems. 

While it seems difficult to simultaneously reproduce the velocity width and column density 
data of damped Lya systems within the context of simple disk models in a CDM cosmogony, our 
model has enough parametric freedom to yield an acceptable fit. In Figure 2 (panel f) and Figure 
3 (the dashed-dotted line) we show such a model with fairly thick disks (k = 0.39), large velocity 
cutoff (u™*^ = 150km s“^), and small neutral fraction (/^r = 0.33). However, such a fit requires 
significant “fine tuning” and yields a physically unmotivated model. As such, we believe that 
simple disk models in a CDM cosmogony cannot reproduce straightforwardly the observational 
data of damped Lya systems at high redshift. This is in contrast to the success of such models in 
explaining the properties of present day spiral galaxies (e.g. Kauffmann 1996). Nevertheless, it is 
interesting to note that CDM models for the formation of structure, combined with assumptions of 
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gas cooling and star formation, tend to significantly overproduce the abundances of faint galaxies 
[White fc Frenk IQQlj . Though there are alternative explanations to this apparent discrepancy 
between observations of the present-day luminosity function and the model predictions, one solution 
would be a suppression of gas cooling and star formation in low circular velocity halos. If cooling 
of gas in less massive Vc < 200km s“^ halos would be similarly suppressed at high redshift, than 
protogalactic disks within a CDM cosmology would remain as a viable model for the nature of 
damped Lya systems at high redshift. 

There are two physical effects which have not been taken into account in our model and which 
could alleviate the discrepancy between the observationally determined velocity width distribution 
and our model predictions. Mo et al. (1997) compute rotation curves for galactic disks within 
realistic dark matter halos. They find that low Vc halos have rotation curves which have a peak of 
approximate height Vrot ~ l.Suc- Since this peak is approximately located at a distance from the 
disk center where damped Lya absorption is dominated (at the threshold for damped systems N^) 
the “effective”rotation velocity of low Vc halos may exceed Vc- Further, hydrodynamic simulations 
employed in the analysis of Ma et al. (1997) indicate that there may be partial ionization of disks 
at larger distances from the disk center. Lines of sight which penetrate the disk at large impact 
parameters typically yield only small Av. If those lines of sight do not contribute to damped Lya 
absorption, the average observed Av of a given disk increases. It remains to be shown if these 
effects can signihcantly improve the comparison between the observationally determined velocity 
width distribution and the disk model predictions. 

In summary, we have compared observationally determined velocity width and column density 
distributions of damped Lya systems at high redshift to theoretically predicted distributions within 
the context of protogalactic disk models for damped Lya systems in a closed CDM model with 
vanishing cosmological constant. Whereas the column density distribution may be reasonably 
well reproduced by assuming rotationally supported disks with scale length comparable to those 
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inferred from successful semianalytic models for the formation of present-day spiral galaxies, 
the simplest protogalactic disk models cannot reproduce the velocity width distribution. This 
discrepancy is mainly due to the large fraction of small circular velocity halos within the CDM 
model. Approximate agreement between the observational data and the disk model predictions 
may be obtained if the contribution of low circular velocity systems Vc < 200 km s“^ to damped 
Lya absorption is significantly diminished. By increasing the thickness of disks and slightly 
decreasing the neutral fraction of gas one may further improve the agreement between observational 
data and the model. Leaving the realm of the model investigated in this paper, there are two 
alternative solutions to the problem. It may well be that only a fraction of all damped systems 
are “regular” protogalactic disks. In this case most of the damped Lya systems which show large 
intrinsic velocity width Av > 100 km s“^ have to be the result of damped Lya absorption resulting 
from lines of sight through double disks, ongoing mergers, and irregular gas distributions. Since 
the observed fraction of damped systems with Av > 100 km s“^ is large (~ 50%) only a similarly 
large fraction of such systems at high redshift z ~ 3 may reproduce the observations. Further light 
on this issue may be shed by hydrodynamical simulations of structure formation [e.g. Haehnelt et 
al. 1997]. Finally, it may simply be that the CDM models employed in this paper overproduce the 
abundances of small Vc halos and do not represent the “correct” structure formation scenario. 

We acknowledge many useful discussions with Martin Haehnelt, Houjun Mo, Simon White, 
and Arthur Wolfe. JXP would especially like to thank W.L.W. Sargent and L. Lu for providing 
their HIRES spectra. This work was supported by NASA grant NAGW-2119 and NSF grant AST 
86-9420443 (JXP) and... 
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4 Figure Captions 


Figure 1 

The doubly conditional cross-section, in units of for a disk yielding a velocity width Av 

greater then a fraction, fc, of the rotation speed for disks with two different central column densities 
{N± = 10^°-^, 10^^-^; top and bottom lines respectively) and thickness (k = 0.108,0.232; dashed 
and solid lines respectively). 

Figure 2 

Optical depth per unit redshift, dr/dz, for damped Ly-a systems showing a velocity width, 
Av, more than a threshold Av > Avth, as a function of Avth- The observationally determined 
distribution is shown by Tier bars for a sample of 26 damped Ly-a systems with mean redshift 
z = 2.87. The dotted lines show theoretically determined velocity width distributions at the 
same redshift for cold, rotating disks in CDM models with Hubble constant h = 0.5 and 
baryonic contribution to the closure density, = 0.05. Panel (a) corresponds to models with 
parameters / = 0.108, k = 0.232, /^r = 1, = 50km s“^, and varying CDM normalization 

as = 2,1,0.6667,0.5 (from top to bottom, respectively), where ag is the linear rms fluctuation 
amplitude on a scale of 8/i“^Mpc, / is the exponential disk scale length in units of the halo virial 
radius, k is the ratio of scale height to scale length, /^r is the fraction of baryons contributing to 
neutral disk gas, and is the lowest circular velocity of halos where gas may still cool and form 
damped Ly-a systems. The remaining panels correspond to models with parameters (b) / = 0.108, 
K = 0.232, as = 1, /at = 1, and varying velocity cutoff = 50,100,200,300 km s“^ (from top 
to bottom, respectively); (c) / = 0.108, us = 1, fjy = 1, u™*” = 50 km s“^, and varying disk 
thickness k = 0.5,0.232,0.108,0.05 (from top to bottom, respectively); (d) k = 0.232, us = 0.5, 
/at = 1, = 50 km s“^, and varying disk size / = 0.834,0.5,0.232,0.108 (from top to bottom 

at high Av, respectively); (e) / = 0.108, k = 0.232, ag = 1, = 50 km s“^, and varying neutral 

gas fraction /tv = 1,0.5,0.25,0.1 (from top to bottom, respectively); and (f) / = 0.108, k = 0.39, 
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as = 1, = 150 km s and /tv = 0.33. 

Figure 3 

Column density distribution for damped Ly-a systems, /(IV), as a function of N at mean 
redshift z = 3. The crosses show observational data taken from Wolfe et al (1995) with the 
vertical lines indicating the Tier ranges. The different lines show theoretically determined column 
density distributions for some of the models in Figure 2. The lines correspond to model parameters 
/ = 0.05, K, = 0.232, cjs = 1, /tv = 1, and = 50 km s“^ (solid line); / = 0.834, k = 0.232, 
as = 2, /tv = 1, and = 50 km s“^ (dotted line); / = 0.108, k = 0.232, as = 1, /tv = 1, and 
^min _ 2 qq g-i (short-dashed line); / = 0.108, k = 0.232, as = 1, /at = 0.1, and = 50 
km s“^ (long-dashed line); and / = 0.108, k = 0.39, as = 1, fjv = 0.33, and = 150 km s“^ 
(dashed-dotted line); 
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